Searching for a Supersolid in Cold Atom Optical Lattices 
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We suggest a technique for the observation of a predicted supersolid phase in extended Bose- 
Hubbard models which are potentiaUy realizable in cold atom optical lattice systems. In particular, 
we discuss important subtleties arising from the existence of the trapping potential which leads to 
an externally imposed (as opposed to spontaneous) breaking of translational invariance. We show, 
by carefully including the trapping potential in our theoretical formalism, that noise correlations 
could prove instrumental in identifying the supersolid and density wave phases. We also find that 
the noise correlation peak width scales inversely with the relative size of trapped Mott domains. 
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Cold atom optical lattice systems are fast becoming 
ideal model systems for studying quantum phases and 
quantum phase transitions in various strong correlation 
condensed matter Hamiltonians 0. Experimental indi- 
cators for novel quantum collective phases such as the 
superfluid and the Mott insulator |^ Q already exist 
in cold atom systems, and there are concrete predictions 
under feasible and realistic experimental conditions, for 
the existence of several other interesting collective phases 
including various quantum spin phases i^] , fractional 
quantum Hall phases 0, and a strongly correlated su- 
persolid cold atom phase 0- 

A key issue in the context of correlated cold atom quan- 
tum phases is their actual experimental observability (in 
contrast to their possible physical existence). The exper- 
imental techniques (based mostly on interferometry) for 
studying cold atom optical lattices differ vastly from the 
standard spectroscopic, thermodynamic, and transport 
measurements prevalent in quantum condensed matter 
systems. Therefore, it is not always manifestly obvious 
how to experimentally observe and compellingly estab- 
lish the various strongly correlated collective phases in 
cold atom optical lattice systems even if their existence 
is convincingly theoretically predicted. The situation is 
further complicated by the existence of the slowly vary- 
ing externally imposed confining trap potential which is 
usually present in cold atom systems. In fact, it is known 
that the confining trap potential introduces special fea- 
tures which, if not accounted for carefully in the theory 
could lead to incorrect conclusions about the quantum 
phase diagram of the system. 

We consider the experimental observability of the pre- 
dicted cold atom supersolid phase u sing the recently pro- 
posed [sl and demonstrated 0, 0, noise correlation 
technique which is particularly suited to probing many- 
body states of cold atoms by directly investigating second 
order correlations in the absorption images of expanding 
atomic gas clouds. Conceptually, the noise correlation 
technique is akin to solid state neutron, electron, or x- 
ray spectroscopy of strongly correlated condensed matter 



systems where a suitable correlation function is spectro- 
scopically observed providing a direct measure of strong 
correlation effects. We show that the noise correlation 
function proves vital in distinguishing between sponta- 
neous and externally broken translational symmetry in 
trapped cold atom optical lattice systems. 

We study the observable properties of two-dimensional 
cold atom, optical lattices in the presence of a parabolic 
trap modeled by the single-band, extended Bose- 
Hubbard Hamiltonian: 
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V ^ UiUj - ^HiUi, (1) 
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where &| creates a boson at the site i and Ui = The 
first term indicates the hopping, t, between nearest neigh- 
bors, denoted by angular brackets. The second term rep- 
resents the on-site interaction energy cost, U. The third 
term contains the interaction energy cost, V, between 
nearest neighbors arising from the spatially extended 
range of interaction possibly generated from a long-range 
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dipolar interaction [l2l Il3| | , a short-range interaction be- 
tween bosons in higher bands of the lattice 0, or a 
combination thereof. A parabolic trapping potential of 
strength k externally breaks translational invariance by 
modifying the chemical potential: fj,i = /xq — 



where Ri — {ix,iy) for a two-dimensional square lattice. 
We work in units of the inter-site spacing, a, equal to half 
the wavelength of the lasers defining the optical lattice. 

To study the experimentally observable ground state 
properties of Eq. |^ we first consider the Gutzwiller, 
mean-field phase diagram in the absence of a trapping 
potential, k = 0]. We obtain the phase digram us- 
ing a variational state 0, in the Fock number basis, 
IV'i) = Y^N^Jo^ /ArJ^i)i where we fix the SN^ vari- 



ational parameters on S sublattices. We find that 
Nm = 6 ensures convergence in the calculations shown 
here. Evaluating mmKtpi^slHlijjies)] yields the phase 



2 




FIG. 1: Mean-field phase diagram of the Bose-Hubbard (ex- 
tended Bose-Hubbard) model in the top (bottom) panel show- 
ing chemical potential versus hopping in an infinite system. 
The Mott and density wave phases are indicated by black 
at integer and half filling, u, respectively. The grey region 
indicates supersolid and the remaining white region super- 
fiuid. The points A through E label the following values of 
{t/U, tJ-o/U): (0.05,0.6), (0.02,0.6), (0.14,0.8), (0.09,0.8), and 
(0.015,0.8). The arrows show the chemical potential shifts 
due to an external trapping potential. 



diagram shown in Fig. ^ In the top panel we choose 
V ^ giving only two phases: the Mott phase and the 
superfluid phase. The incompressible Mott phase occurs 
at integer filling, i/, with order parameter (rii) Sinteger 
while the superfluid can have variable filling with or- 
der parameter (6^) 0. The bottom panel shows the 
mean-field phase diagram with V = U/z, where 2 = 4 
is the number of nearest neighbors. Two new phases ap- 
pear: the incompressible density wave phase at = 1/2 
and the supersolid phase. The density wave phase on 
a square lattice has checkerboard oscillations in density 
with order parameter (— 1)(*^+'!'^ ((n^) — v). The super- 
solid consists of coexisting density wave and superfluid 
order. Monte Carlo studies suggest that the delicate su- 
persolid phase successfully competes with phase separa- 
tion 16] for y > J7/z. We seek distinct experimental in- 
dicators of density wave and supersolid order as captured 
by a zero temperature, modified Gutzwiller mean-field 
theory in a model which includes the trapping potential. 
The technique employed here can be extended to include 
finite temperature effects. In this work we focus on zero 
temperature properties as a starting point. 

Inclusion of the trapping potential mixes phases within 
the trap and may obscure experimental signatures of the 
phases discussed above. The center of the trap corre- 



sponds to a point in the phase diagrams of Fig. ^ The 
arrows indicate positions moving towards the edge of the 
trap which see an effectively lower chemical potential 
thereby mixing phases within the trap. We incorporate 
mixing by modifying the Gutzwiller ansatz to bring out 
an important realistic feature; smooth modulation of the 
density throughout the trap. We allow the variational 
parameters to vary from site to site while minimiz- 
ing the energy of the total Hamiltonian: min [(^Pli? j^*)] , 
where j^") = Yii IV'i) generalizes jf/'i)- This leaves a mini- 
mization problem over NsNm variables, where Ns is the 
total number of sites. As an initial configuration we 
choose the output from our initial mean-field minimiza- 
tion performed independently at each site. First order 
correlation functions calculated with this two step proce- 
dure show ^3 close agreement with Monte Carlo studies 
[lij for = and away from the phase transition. This 
suggests that we retain quantitative accuracy for 7^ 
although the techniques employed here remain untested 
against small system Monte Carlo results for V 0. In 
what follows we consider system sizes comparable to ex- 
periment Ns = 101 X 101 with k/U = 7 X 10"'' yielding 
a total number of particles TV ^ 10^ for hq/U — 0.6 — 0.8 
and7V = ^,(*|n,|*). 

Time of fiight absorption imaging of atoms released 
from trapped optical lattices provides a direct probe 
of first and second order correlation functions of the 
initially trapped state. We first consider the momen- 
tum distribution function measured through an average 
of images of expanding atomic clouds after a time T. 
Assuming that the particles interact weakly after the 
trap is adiabatically turned off, the number of particles 
at position r in the expanding cloud is given by |l8l |: 
{n{r))T ~ \w{Q{r))\'^ p{Q{r)), where Q{r) = amr/HT 
for particles of mass m. The angular brackets indi- 
cate averaging, which, experimentally, occurs over sev- 
eral shots, w is the Fourier transform of the Wannier 
function. For illustration we take particles in the low- 
est band of a deep lattice. The harmonic ^proximation 
then gives (using the conventions of Ref. 8]): [^(Q)^ « 
|w(0)|2exp[-|Q|7(7r2/Vi)], for a lattice depth Vl in 
units of the photon recoil energy, p is a first order corre- 
lation function related to the distribution of lattice mo- 
menta k in the trapped system: 

p(fc) = ^e'(«-«^)-'=(6t5^.). (2) 

Peaks in p indicate inter-site phase coherence. From shot- 
to-shot a coherent state maintains, to a good approxima- 
tion, the same phase relationship between Fourier compo- 
nents in the averaging process. However, for incoherent 
states with random Fourier components, i.e. the Mott 
and density wave phases, the averaging process elimi- 
nates any peak in the Fourier spectrum. One can over- 
come this loss of information by averaging the second 
order, noise correlation function. 
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Noise in the imaged atomic momentum distribu- 
tion provides a direct measure of higher order corre- 
lation functions. Shot-to-shot averaging in the quan- 
tity (n(r)n(r'))r — (n(r))r(n(r'))T yields results pro- 
portional to momentum correlations in the ground state 
of the trapped system 0: 

g{r,r') ^ \w{Q{r))\^\w{Q{r'))\^G{Q{r),Q{r')) 

and, in contrast to the first order correlation function, 
does not vanish in finite sized system measurements of 
incoherent states. Higher order terms arise in the nor- 
malized, second order correlation function, G{k,k'), de- 
fined as: 



G{k,k') 
p{k)p{k') 



y e*(«».' 6^5. ,6-,) 



p{k)p{k') 

ma\ 5{k — k') 
W) |w(fc')pp(fc')' 



- 1 
(3) 



providing a direct measure of Mott and density wave or- 
der. The delta function results from normal ordering 0. 
We evaluate p and G using the output of our modified 
Gutzwiller mean-field theory, 

In a translationally invariant system (i.e. k = 0) p and 
G qualitatively distinguish the ground states of Eq. Q . 
Macroscopic occupation of a Bloch mode of the optical 
lattice at the reciprocal lattice vectors should lead to 
peaks in p demonstrating a spontaneously broken U{1) 
symmetry associated with the global phase, and there- 
fore superfluid order. If the peak exists at wavevectors 
other than the reciprocal lattice vectors we expect, in 
addition, spontaneously broken translational symmetry, 
i.e. supersolid order. (In contrast, applied potentials 
with spatial and temporal j^^l modulation of the 
the host lattice have led to the observation of coherence 
peaks away from the reciprocal lattice vectors not asso- 
ciated with a supersolid.) Furthermore, strong peaks in 
G at the reciprocal lattice vectors of a translationally in- 
variant optical lattice indicate Mott order while strong 
peaks in G at half the reciprocal lattice vector indicate 
spontaneously broken translational symmetry and there- 
fore density wave order. As we will see both correlation 
functions, in a trapped system, merely quantitatively dis- 
tinguish between the supersolid and other phases. 

We first verify that, as expected, superfluid and 
Mott order in the Bose-Hubbard model remain identi- 
fiable through first and second order correlation func- 
tions in a trapped system. Fig. |21 plots both correla- 
tion functions for k/U = 7 x 10"'' at parameters cor- 
responding to the points A (top panel) and B (bot- 
tom panel) of Fig. ^ as a function of k. The left y- 
axis plots the normalized momentum distribution func- 
tion, p{k) = p(fc)|?i)(fc)p/|iVi5(0)p, while the right y-axis 
plots G{k,0). The superfluid phase (top panel) shows a 
strong coherence peak at the reciprocal lattice vectors 
K = (27ri,2^j), for integer j, (p(27r, 2^)/p(0, 0) - 0.3) 




FIG. 2: The first order, p, (left y-axis, dashed line) and 
second order, G, (right y-axis, solid line) normalized correla- 
tion functions plotted as a function of wavevector, k, for the 
trapped Bose-Hubbard model. We set k = kx = ky. The 
dashed lines are shifted up by 0.02 and right by 7r/4 for clar- 
ity. The superfluid phase (top panel, point A in Fig.0 shows 
strong peaks in p at the reciprocal lattice vectors K . In this 
and in the following flgures we choose a deep optical lattice, 
Vl = 20, which determines the weight of p at higher k. The 
Mott phase (bottom panel, point B in Fig. shows strong 
peaks in G at the reciprocal lattice vectors K . The inset plots 
the full width at half maximum of the central peak in G at 
k — (27r, 27r) versus hopping for fi/U = 0.7. The vertical ar- 
row indicates the position of a superfluid to Mott transition 
at the trap center. 



with little signal in the noise. The exponential decay of 
w(fc) strongly suppresses the peak in p at {2tt, 2tt) as op- 
posed to the peak at {2tt, 0) (not shown). The Mott phase 
(bottom panel) is identifiable through a strong peak in 
G, at the reciprocal lattice vectors. Note that here there 
is still a small amount of superfiuid at the edge of the 
system, identifiable through weak peaks in p. 

Individual noise correlation peaks reveal position in 
the phase diagram through an additional feature, peak 
width. Increasing t/U through the Mott phase we find 
the noise correlation peaks to decrease in height and split. 
The splitting results from our choice of normalization. 
Other suitably chosen normalization schemes can elimi- 
nate the splitting by excluding coherence peaks from the 
definition of the normalization |2H . In addition, we pre- 
dict that the widths of the peaks vary inversely with the 
size of the Mott domains (rather than the total system 
size), in qualitative agreement with recent experimental 
results j21|. The inset to Fig. El shows the fufl width 
at half maximum along the kx axis, Akx, of the central 
(27r, 27r) noise correlation peak increasing with t/U. An 
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FIG. 3: The same as Fig.|21but for the points: C (top panel, 
superfluid), D (middle panel, supersolid), E (bottom panel, 
density wave) of Fig. 



abrupt increase in peak width occurs when the center of 
the trap crosses from the Mott to the superfluid regime. 
We conclude that noise correlation peak widths yield de- 
tailed information related to the location of the system 
in the phase diagram. 

Fig. 121 shows correlation functions of the extended 
Bose-Hubbard model {V/U = 0.25) along the points C 
(top panel), D (middle panel), and E (bottom panel) of 
Fig-ffl The top panel shows signatures of superfluidity; a 
strong peak in p, with p{2tt, 27r)/p(0, 0) ~ 0.18, and weak 
structure in G. The bottom panel, deep in the density 
wave regime, shows signatures of checkerboard density 
wave order only; strong peaks in noise correlations at 
k = K and K /2. The paucity of superfluid at the edge 
yields almost no structure in p. The middle panel is plot- 
ted for parameters with supersolid at the trap center and 
superfluid near the edge. Here we see a coherence peak 
at the reciprocal lattice vectors K and very weak struc- 
ture in G. There is also a very weak peak (~ 0.01) in p 
at (tt, tt). We conclude that here weak supersolid order 
is nearly indistinguishable from superfluid order. 

The supersolid shows a small signal in our mean-field 
estimates of the first order correlation function. We 
define the strength of a supersolid phenomenologically 
through p{'n, tt)/p{0, 0) which is sensitive to spatial mod- 
ulation of the superfluid order parameter, (bi). The 
nearby superfluid and density wave phases surrounding 
the supersolid phase in parameter space show no modula- 
tion in {bi). We then expect the supersolid to exhibit only 
weak structure in (bi) and therefore p(7r,7r). The mod- 
ulation amplitude can be enhanced if the density wave- 
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FIG. 4: The main panels plot the same as Fig. |5| but for 
V/U — 1 and fi/U = 0.8 with the dashed lines shifted up by 
0.01 and right by 7r/4 for clarity. The main top (bottom) panel 
is in the supersolid (density wave) regime with t/U = 0.2 
{t/U = 0.1). The insets show the phase diagram (as in Fig.Q 
but with V/U — 1) near the 1^ — 1/2 density wave. 



supersolid and superfluid-supersolid phase boundaries lie 
far apart in parameter space. We may therefore reach the 
strong supersolid regime by increasing V/U. 

We now consider the strong supersolid regime in a trap. 
The main panels of Fig. 01 plot the same as Fig. [3 but 
for V/U = 1, fi/U = 0.8, t/U = 0.2 (top panel), and 
t/U — 0.1 (bottom panel). We choose the parameters 
for the top panel so that the supersolid is strongest at 
the center of the trap with superfluid near the edge as 
indicated in the inset. We again see a strong, first order 
coherence peak a.t k — K with weak structure in G. We 
also find p(7r, 7r)/p(0, 0) - 0.03. 

We now compare the correlation functions of a strong 
supersolid with those of the trapped density wave. The 
bottom panel of Fig. ^ lies in the density wave regime 
with a significant amount of superfluid near the edge. A 
strong enhancement in noise correlations at K /2 signals 
the onset of density wave order. Superfluidity near the 
edge enhances p{K) and splits the G{K, 0) peaks. There 
is also a small but surprising feature in 'p near (tTjTt). 
The weak (tt, tt) feature, p(7r, 7r)/p(0, 0) ~ 0.02, previ- 
ously associated with the supersolid phase, arises here 
as the density wave and trap corrugate the superfluid 
at the edge. The modulated superfluid does not spon- 
taneously break translational symmetry along the edge 
but consists of a mixture of phase boundaries mitigated 
by the trapping potential. We further find that the edge 
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mode density oscillations never reach zero, as in the case 
of a bulk supersolid. The essential differences in going 
from the top panel (comprising a trapped supersolid) to 
the bottom panel (a trapped density wave) are therefore 
merely quantitative; a large enhancement in G{K /2,0) 
combined with a sharp drop in p{K/2) and p{K). 

Our interesting and important findings are that (1) 
noise correlations enable relatively precise determination 
of phase boundaries in trapped systems and are there- 
fore fundamental in providing unambiguous observation 
of the supersolid phase, and (2) the edge of the cold atom 
density wave may accommodate a superfluid component, 
mimicking the existence of an edge supersolid. 
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